Polarons in 2+e dimensions and giant figure of merit of semiconducting nanolayers 
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Polarons - electrons coupled with lattice vibrations - play a key role in the transport and optical 
properties of many systems of reduced dimension and dimensionality. Their confinement affects 
drastically the phonon, polaron, bipolaron and multi-polaron properties of quantum wells. Here 
we calculate the energy spectrum and thermopower of polarons confined to a potential well as 
a function of the well thickness. We show that the polaron mass enhancement in 2+e dimensions 
explains a giant thermolectric power recently observed in doped semiconducting nanolayers (multiple 
quantum wells (MQWs)), and propose a route for enhancing the performance of thermoelectric 
energy nanoconverters increasing their figure of merit by more than one order of magnitude. 



o 
o 



> 

(N 

(N 
&\ 

o 



X 



While the form of the electron-phonon (Frohlich) 
Hamiltonian [l[ is the same in all dimensions, the po- 
laron mass, binding energy and different response func- 
tions depend on the dimensionality for large polarons [3[ 
and on the dimensionality and crystal lattice structure 
for small polarons Q , especially in the intermediate cou- 
pling regime. In particular, directly following |5J], the 
variational large polaron energy was calculated in Ref. Q 
in different dimensions providing scaling relations valid 
to the second order in the Frohlich coupling constant a. 
In the important case of a two-dimensional confinement 
(2D), scaling takes the form Q: 
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for the ground state binding energy (polaron level shift, 
E p ), polaron mass, to*, and the mobility, /i (if the scat- 
tering is dominated by optical phonons). Here m is the 
band mass in a rigid lattice. 

These scaling relations are instrumental for a qualita- 
tive understanding of physical properties of semiconduct- 
ing layered nanostructures with polaronic carriers, which 
have a potential to provide a new route for realizing prac- 
tical thermoelectric devices with the high figure of merit, 
Z = S 2 ct/k (S, a and k are the thermoelectric Seebeck 
coefficient, and the electrical and thermal conductivities, 
respectively). Following a theoretical prediction by Hicks 
and Dresselhaus @ that layering can increase the figure 
of merit due to an increase of the density of states of con- 
fined carriers, different MQWs have been manufactured 
with an enhanced Z. More recently Ohta et al. [3] ob- 
served a giant enhancement of S, in Nb:SrTi03/SrTi03 
superlattices consisting of conducting Nb-doped SrTiC>3 
(STO) nano-layers confined between insulating STO lay- 
ers. With decreasing thickness of the conducting layer, 



the magnitude of S increases to a value about five times 
that of the bulk Nb:STO, reaching almost 500 fiV/K for 
the one unit cell layer, which is much higher than what 
was anticipated theoretically @. 

Doped STO has been recognized as an archetypical po- 
laronic perovskite a long time ago Q. Some earlier and 
recent experimental observations of the optical conduc- 
tivity in the bulk Nb-doped STO 0] clearly reveal the 
evidence of the mid-infrared optical conductivity band 
provided by the polaron mechanism as in many other ox- 
ides. The effective mass of the charge carriers is obtained 
by analysing the Drude spectral weight. Defining the 
mass renormalization of the charge carriers as the ratio of 
the total electronic spectral weight and the Drude spec- 
tral weight, a two-fold mass enhancement is obtained, 
which is attributed to the electron-phonon interaction 
(EPI). The missing spectral weight is recovered accord- 
ing to the sum rule [lfj in the mid-infrared optical con- 
ductivity band. This band results from EPI, tradition- 
ally associated with the polaronic nature of the charge 
carriers. The effective mass obtained from the optical 
spectral weights yields an intermediate electron-phonon 
coupling strength, a ~ 3 @. Importantly, very recently 
Choi et al. ll| using optical spectroscopy observed a 
three-fold mass enhancement in Nb:SrTi03/SrTi03 su- 
perlattices compared with the bulk Nb-doped STO as the 
Nb:SrTi03 layer thickness decreased from eleven to one 
unit cell. The authors of Ref. [ll| have suggested that 
such an increase should be attributed to a dimensional 
crossover of a polaron from 3D to 2D, and seem to be as- 
sociated with the strongly enhanced Seebeck coefficient. 
Here we propose a theory of the dimensional crossover 
of the large polaron from 3D to 2D providing detailed 
explanation of the giant thermopower enhancement in 
polaronic MQWs and predicting further possibilities for 
enhancing the performance of thermoelectric MQW de- 
vices. 

We consider a polaron confined to an infinite potential 
well of the width d in z-direction, such as Nb:SrTi03 
layer merged into a bulk undoped SrTiOs. For large (or 



2 



continuous, in other words) polarons the lattice structure 
is irrelevant (2|, so that we adopt an effective band mass 
approximation and the canonical Frohlich interaction [l| 
of doped carriers with optical phonons in ionic insulators, 
and also include the confinement potential U(z): 
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Here d q (d^) annihilates (creates) optical phonons with 
the frequency to. While EPI is moderate, one can apply 
the perturbation theory with respect to Vq, which works 
well even in the intermediate coupling regime a < 2 as 
verified by the Feynman all-coupling path-integral the- 
ory H ■ The quantum states of a noninteracting electron 
and phonons are classified with respectively the in-plane 
2D electron wave- vector, kn , and a subband index n, and 
the 3D phonon wave-vector q. The unperturbed state is 
the vacuum |0) of phonons and the electron plane wave 
|kii,n, 0) = -^e lk H' P |n)|0), where 7? is the in-plane elec- 
tron coordinate, and \n) depends on z. The bare en- 
ergy spectrum of the electron in a layer of the area L 2 

h 2 k 2 

is £7 n k|| = ~2^r + e n, and n — 1,2,..., corresponds to 
quantized subbands due to confinement in z direction, 
perpendicular to the layer. EPI couples the state ) with 
the energy E n ^ and states of a single phonon with the 
momentum q and the electron with the energy E n rk»—q« ■ 
The corresponding matrix element is 



(k|| — qt|| , 77.'; lq| Vq|k|| , ?7.; 0) = 
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where f2 is the volume of the whole sample including the 
doped layer and the host bulk crystal. The renormalized 
second-order energy E n ^,, is found as 
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We consider temperatures ksT small compared with 
the characteristic phonon energy fkv and low densities 
of carries. In this case there is no imaginary part of 
-E^iky for sufficiently low energy subbands n, which means 
that the momentum is conserved. Expanding over k\\ and 
replacing the sum over n' for an integral, one obtains 
£"k = —ahiu + h 2 k 2 /2m3£> in the 3D system, d — * oo, 

m 3D = m/(l- 
n' = 1 and the 
0), we find E k — 



with the familiar 3D polaron mass 
a/6). Keeping only terms with n 
matrix element 



s *9»*| n ) = 1 in Eq.i 
—TTaHu/2 + h 2 k 2 /2rri2D with the polaron mass m* 2B = 
m/(l — an/ 8) for a 2D layer, d — > 0. Both the polaron 
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FIG. 1: Normalized polaron level shift E pn /ahuj (top panel) 
and the relative polaron mass renormalization (1 — m/m*)/a 
(bottom panel) for a number of subbands n and different 
nanolayer thickness t in the crossover region. A sudden 
increase in the shift and in the mass for sufficiently thick 
nanolayers ( t > 2) and high energy subbands is due to a 
breakdown of the perturbation expansion. 



binding energy and the mass satisfy the scaling relations 
Eqs.©. 

In the crossover region, where neither 3D nor 2D limits 
are applied, the spectrum is quantized into subbands, 
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with the subband polaron energy shift —E pn and the sub- 
band in-plane polaron mass, to* , found as 
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respectively. Here, we use dimensionless phonon mo- 
menta y — q z d/n, x — q\\d/n, z-coordinate u = nz/d of 
the polaron confined to the layer of the thickness d, and 
dimensionless thickness of the layer, t = q p d/n. When 
polarons are confined to the infinite potential well within 
< z < d, the bare eigenfunctions and bare eigenvalues 
are |n) = (2/d) 1 ^ 2 sm(irnz/d) and e„ = h 2 ir 2 n 2 / (2md 2 ), 
respectively. Then, integrating first over the transverse 
momentum y to avoid complicated oscillating functions 
in Eqs.dHE]) yields 
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for the off-diagonal terms, n' ^ n, where m = n' + n and 
k = n' — n (n',n = 1,2, ...oo). Two last terms in the 
diagonal contribution, Eq. (fi"2")) , and the last term in the 
off-diagonal contribution, Eq. (ll3|) are numerically small 
about a few percent of the remaining terms. Integration 
over x and summation over n', performed using any stan- 
dard software, yield the shift and the mass interpolating 
well between 2D and 3D limits. However, the polaron 
level shift and the polaron mass are different for different 
subbands in 2 + e crossover region as shown in Fig[TJ As 
can be seen in FigtH one encounters the familiar flatten- 
ing of the polaron dispersion and the breakdown of the 
perturbation theory at energies about the phonon energy 
[l2j for high-energy subbands. 

Importantly, the polaron mass changes from its 2D 
value to the 3D value for t in the range of about unity, 
Figffl The experimental values of the effective thick- 
ness t are about t — 0.16JV, where N is the number of 
atomic layers in the well 0, [ll|, if we take fko = 0.05 



eV, m = m e , and the lattice constant a — 0.4 nm, which 
are typical of STO. Hence, the experimental conditions 
are just right for a significant change of the polaron mass 
with the number of layers in Nb:SrTi03/SrTi03 super- 
lattices where N changes from 1 to 11 @, [HI]. The 
three-fold mass enhancement with decreasing thickness 
[Tl| corresponds to a w 2.8, FigJ^l which is close to opti- 
cally measured value in the bulk Nb:SrTiO @. 

Equations (jTOj) and (fTTjl and the Boltzmann theory al- 
low us to analyze the dependence of different kinetic and 
thermodynamic properties on the confinement. In par- 
ticular, applying the energy- independent relaxation time 
(t) approximation, one obtains the Seebeck coefficient 
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where f{E) = 1/[1 + exp(£' - Q/k B T] is the Fermi- 
Dirac distribution function. In the bulk system with 
m^ D = 2m e , the Fermi energy is about C,/k B = 1250K for 
the polaron density n p = 10 21 cm -3 , so that carriers are 
almost degenerate at room temperature. However, the 
system evolves from a degenerate to a classical heavier 
polaron gas with decreasing nanolayer thickness because 
of the strong mass enhancement in the crossover region, 
Fig. 3. Hence one has to use the exact Fermi-Dirac statis- 
tics in this region. Integrating over momentum in Eq. (fT4|) 
yields: 
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where y n = exp ^ - ^^J- 
following sum rule 
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satisfies the 



(16) 



and Li2(z) = J2^ =1 z k /k 2 is the dilogarithm function. 
Since Huj/k B T > 1 at room temperature for typical opti- 
cal phonons in oxide ionic lattices, one can safely ne- 
glect all higher subband contributions in Eqs. lfTSl [T6|) 
compared with the lowest subband term n = 1 in the 
crossover region, t < 1, avoiding the breakdown of the 
perturbation expansion, FigfT]. The thermopower, calcu- 
lated using this single-subband approximation, accounts 
rather well for the experiment pj with the same values 
of the EPI coupling, which accounts for the experimen- 
tal polaron mass in FigJ^l The carrier density turns out 
significantly lower compared with the chemical density 
presumably due to a partial localization of polarons by 
impurities. 

Finally, let us discuss the confinement enhancement 
of the figure of merit. In semiconducting MQWs such 
as Nb:SrTi03/SrTi03 superlattices, the thermal conduc- 
tivity is mainly due to phonons and largely indepen- 
dent on doping and confinement. If we take S oc m* 
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FIG. 2: (a) Polaron mass m* /m of the lowest subband 
(top panel) and Seebeck coefficient S (bottom panel) in the 
crossover region for a few coupling constants a and hu = 0.05 
meV, T = 300K, m = m e and n v = 1.66 x 10 20 cm" 3 (solid 
lines) as the functions of the nanolayer thickness d (inset in 
the top panel). Symbols correspond to the experimental mass 
measured optically in Ref. and to the experimental See- 
beck coefficient measured in Ref. [3]. 



then, using the scaling Eqs.|T|), the maximum confine- 
ment enhancement of Z is estimated as Z 2 d jZ^D = 
TO 3i)(37ra/4) 2 /i3o(37ra/4)/TO3 D (a) 2 ^3_D(a), if polarons 
are scattered by absorption of optical phonons. In this 
case, the mobility depends on the coupling strength as 
[H[ M3£> oc [l — a/6 + 0(a 2 )] /a, so that the enhance- 
ment, Z 2D /Z 3D = (4/37r)(l-a/6)/(l-7ra/8), is only 1.3 
for a = 2. However, at temperatures well below the char- 
acteristic optical phonon temperature hw/kB, the num- 
ber of thermal phonons is exponentially small, so that the 
polarons are mainly scattered off neutral impurities since 
the impurity potential is screened by heavy polarons at 
substantial doping. The scattering time off neutral im- 
purities increases with the polaron mass as r oc (m*) 2 
so that the mobility is proportional to the mass, 
fi oc m*. In this case the confinement enhancement of 
the figure of merit may be giant. For the coupling a = 2, 
Z 2 d/Z 3D = (1 - «/6) 3 /(l - 7ra/8) 3 « 30 is much larger 
than in a free electron layered system Q . This prediction 
provides a further route for enhancing the performance 
of thermoelectric energy converters. 

To conclude, we have developed the theory of po- 
laron crossover from 3D to 2D. In the crossover re- 
gion the polaron energy spectrum is quantized into size 
subbands with the subband binding energy and po- 
laron mass strongly dependent on the size of the con- 
finement nanolayer, when the thickness of the layer is 
comparable or less than the large polaron radius, d < 
l/q p . We have shown that this condition is satisfied in 
Nb:SrTi03/SrTiC>3 superlattices and proposed a detailed 
explanation of the giant thermopower enhancement ob- 
served in these structures 0, [H| • The theory predicts a 
giant confinement enhancement of the figure of merit in 
those polaronic MQWs where the scattering of polarons 
is dominated by impurities. 
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